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Abstract 

The higher twist effects in deep inelastic scattering on the nuclear target are 
studied using time ordered perturbation theory. We showed that the collinear rescat- 
tering of the outgoing quark on the extra nucleons via the contacting gluon-pair is 
dominant in nuclear size-dependent effects. The Qiu-Vitev resummation is proved 
by using the geometric properties of the higher twist amplitudes. The leading contri- 
butions of nuclear-enhanced effect to the DGLAP evolution equation are resummed 
in the same framework. 



PACS numbers: 13.60.Hb; 12.38.Bx. 

keywords: Higher twist effect; PQCD evolution equation 



1 



1 Introduction 

Deep inelastic scattering (DIS) experiments are the main source extracting the infor- 
mation about QCD structure of the hadronic target. The DIS amplitudes can be expanded 
according to their twists. The leading-twist amplitudes describe the quark-gluon struc- 
ture of a target at short distance, while the higher twist effects and their summation may 
provide a valuable information for the multi-parton correlations at long distance in the 
target [1]. The nuclear target is an ideal "laboratory" for testing higher twist theory 
because of its special multi-parton configuration. 

However, it is difficult to sum all contributions in the twist expansion due to the 
complicated structure of the higher twist amplitudes and the unknown multi-parton cor- 
relation functions. Instead of a complete summation, the resummation of the leading 
contributions to a special effect from the twist expansion is useful. For example, the 
Glauber model is used to resum the contributions of the multi-parton scattering to the 
shadowing effect in the nuclear structure functions in the Glauber- Mueller formula [2]. 
This formula assumes that the leading contributions are from the incoherent scattering 
of color dipole on the nuclear gluons. Another example in this realm is the Jalilian- 
Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation [3], where the 
soft nuclear gluons are treated as a classical color potential because of its large density 
at small Xb- In this background, energetic partons in the multi-dipole configuration in- 
teract with soft gluons through the Wilson lines, which resum the contributions of the 
soft nuclear gluons at the eikonal approximation [4,5]. Both of the above two researches 
are broadly used to predict the saturation phenomena in the nuclear structure functions, 
where unitarity is restored. 
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On the other hand, a different resummation of the nuclear higher twist amplitudes is 
recently proposed by Qiu and Vitev (Qiu-Vitev resummation) [6]. This model assumes 
that the leading contributions to the nuclear size effect in the structure functions originate 
from the singularities of the gluon correlation functions in the higher twist amplitudes [7] . 

The higher twist amplitudes have complicated color-spin structure. The resummation 
of a higher twist effect is valid only if the contributions from a lot of irrelevant diagrams 
can be entirely neglected. Unfortunately, this point is unclear in the researches mentioned 
above. 

The motivation of the present work is to study the geometric properties of the higher 
twist amplitudes for the DIS on nuclear target using time ordered perturbation theory 
(TOPT) [8], which was developed in the collinear approximation (CLTOPT) in [9]. Com- 
paring with the standard covariant perturbation theory (CVPT), the CLTOPT provides a 
more intuitive geometric character of the higher twist amplitudes. Using the above TOPT 
method we study the DIS processes on the nuclear target at arbitrary twist. We present 
a geometric higher twist effect in the DIS on the nuclear target: the leading contribu- 
tions to a nuclear size- dependent effect in the DIS processes originate from the multiple 
rescattering of the outgoing quark on the extra nucleon via the exchange of the contacting 
gluon-pair. It reproduces the Qiu-Vitev resummation [6] but using different dynamics. 
In particular, this space-time character of the amplitudes in the CLTOPT allows us to 
distinguish the leading diagrams from a lot of sub-leading diagrams without calculations 
of the Feynman diagrams. We resum the corrections of the geometric higher twist effect 
to the structure function and DGLAP evolution equation [10,11]. 

The paper is organized as follows. In Section 2 we present a short review about 
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TOPT in the collinear approximation, which is necessary for understanding this work. 
For a complete calculation of the higher twist amplitude including the matrix element, 
we discuss in Section 3 the normalization of the correlation functions. In Section 4, 
using CLTOPT we demonstrate the space-time picture of the higher twist amplitudes, 
which exposes obvious nuclear size-dependence. In Section 5 we derive the expressions 
for the nuclear dependence of the structure function at twist-4 level, then the results are 
generalized to sum the contributions of all twist amplitudes. In Section 6 we study the 
corrections of the above mentioned higher twist effect to the DGLAP evolution equation. 
The discussions and a summary are presented in Section 7. 
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2 Definitions and tools 

We consider DIS of leptons from an unpolarized nuclear target. We denote the light- 
like vectors along the positive and negative directions as = l/\/2(l, 0j_, —1) and n M = 
l/y/2(l, 0j_, 1). Therefore, we have 

n-n = n-n = 0, n-n — 1. (1) 

As is well known, the infinite momentum frame is equivalent to the light-cone coordi- 
nate. Here, however, we would rather use TOPT in the infinite momentum frame than 
the light-cone quantization. The reason is that the contributions of the backward (instan- 
taneous) propagator in the latter case are neglected and it is not a complete description 
for the higher twist processes. In the concrete, we choose an infinite frame, in which the 
target momentum and virtual photon momentum have the following forms, 

Q 2 

q» = q + n + q n = -x B P + n^ + n M , (2) 

ZX B r 

where Q 2 = —q 2 and x B is the Bjorken variable. 

In this work we take the physical axial gauge: n • A — 0, with A being the gluon field. 
The inclusive cross-section for the unpolarized DIS is determined by the hadronic tensor 
W^ u , which can be decomposed into the transverse structure function F t (xb,Q 2 ) and 
longitudinal one F l (xb,Q 2 ), 

W^(x B , Q 2 ) = e L v F L (x B , Q 2 ) + et?F T (x B , Q 2 ). (3) 
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The projection operators are defined as 



< - ^ + S-^) ( xbP+ * + S* n j • (4) 

= hjfn v + rfn» - g^) = ^dT '. (5) 



Clearly, using 



^ • = 1 ef • ef = 1 and ef • ef = 0, (6) 



one can directly extract F t {xb,Q 2 ) and F l (xb,Q 2 ) from according to 



Ft(xb, Q 2 ) = 2e>?W^(x B ,Q 2 ), (7) 
Fl(z b ,Q 2 ) = 4efW^(x B ,Q 2 ). (8) 

Note that the hadronic tensor W 1X111 {xb, Q 2 ) can also be decomposed into the structure 
functions F 1 (x B ,Q 2 ) and F 2 (x B ,Q 2 ) and their relations with Eqs. (7) and (8) are 

F L (x B ,Q 2 ) = F2{XB ' Q2) -2F 1 (x B ,Q 2 ), 
x B 

F T (x B ,Q 2 ) = 2F 1 (x B ,Q 2 ). (9) 

To further examine the role of higher twist corrections in the DIS processes, we present 
the TOPT-form of the QCD-propagators in the collinear limit [9]. As will be discussed 
in details below, the nuclear dependence on both the structure functions and evolution 
equations is dominated by the higher twist amplitudes. A straightforward reason is that 
the amplitude at higher twist contains the multi-parton operators in the matrix element, 
which can come from the different nucleons of a nucleus and give rise to the nuclear 
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size-dependence. However, we have a lot of amplitudes corresponding to the multiple 
interactions among the partons. The important question is: which of them contributes 
the leading nuclear effect to the processes at a given order of approximation? We find 
that the TOPT in the collinear limit, which is developed in our previous works [9], is a 
useful tool for answering the above question. 

TOPT is exactly equivalent to the method of the Feynman diagram [1]. A covari- 
ant propagator in the Feynman diagram is decomposed to the forward- and backward- 
propagators in TOPT and this leads to two remarkable properties: the interaction is 
ordered on time and the propagating momentum is on the mass shell. As we will point 
out in the following, the former shows how a process evolves in the light-cone coordinates 
and gives us information about the nuclear effects, while the latter is useful to decompose 
a complicated calculation. 

At first sight it seems like that TOPT leads to a proliferation of diagrams. However, 
we find that there is only one TOPT diagram corresponding to one Feynman diagram in 
the collinear approximation, where the transverse momentum of the parton is neglected. 
In fact, the TOPT-propagators in the collinear limit have the following special forms [9] 

s F (k) = do) 

for quarks, and 



Gf(k) = 
Gf(k) = 



df 

k ■ n(k • n — k -n) 

n a n f3 

{k-n) 2 ' 



(12) 
(13) 
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for gluons, where "F" and U B" label the forward- and backward-moving partons, respec- 
tively. Using helicity conservation and Eq. (1), the non- vanishing vertices of QCD and 
QED in the collinear approximation are shown in Figs. 1 and 2, respectively. According 
to those selection rules, we draw the twist-4 amplitudes, which contribute to Ft(xb,Q 2 ) 
and F l (xb, Q 2 ) in the collinear limit, in Figs. 3 and 4, where we label the outgoing quark 
or outgoing gluon as "5" or "F", since it contributes ~ 7 • n or ~ d°f, respectively. We 
did not draw the conjugate diagrams. One finds that only one TOPT graph corresponds 
to one covariant Feynman diagram in the collinear limit. 

The deep inelastic interaction in the TOPT-form is time-ordered in a given infinite 
momentum frame and shows how a system evolves in the light-cone coordinates. For 
illustrating this picture, we take the Fourier transformation of the forward- and backward- 
propagators. Note that the right hand sides of Eqs. (10) and (11) indicate that Sf is 
only a function of k~, while Sb is related to k + . Thus, we have 



fF{yt-yi)Kyi-V2)^\yi-yi) 
s F (vt-vt), (14) 
Mvr - V2Myi - yt)^\y{ - vi) 
S B {yi-y 2 )- (15) 

These expressions show that the forward- and backward-propagators for quarks in the 
collinear approximation are parallel to y + - and y~-axes, respectively. Using this property, 
for example, we draw the twist-4 amplitudes of Fig. 3a in the y + — ^"-coordinate in Fig. 
5a. One finds that Sf and Sb correspond to the contact- and instantaneous-interactions 
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J S F (k )exp[ik- (yi - y 2 )\ 



d A k 
(2k? 



J S B (k + ) exp[ik ■ {y x - y 2 )\ 



d A k 
(2k? 



in the light-cone coordinates. The backward propagator of the gluon (Eq. (13)) has the 
same property as Sb, i-e., Gb(Ui — y%) is always parallel to the y~-axis. We focus on 
the relation between the interactions of 7*-parton and the structure of the correlation 
functions, which are defined on the |/~-axis in the collinear factorization. Therefore, for 
simplicity, we project the interactions of 7*-parton in Fig. 5a on the y~-axis (Fig. 5b). 
Thus, we get the geometric diagrams Figs. 6 and 7, corresponding to Figs, 3 and 4, 
respectively. Here we note that Eq. (12) indicates that G F (k) is a function of both k~ 
and k + . Therefore, the forward propagator of gluon has two components — along y + and 
y~ — even in the collinear limit. An exception is that the gluon propagator connects with 
the hard- and soft-parts of the amplitude. The factorization scheme [12] considers that 
soft gluons attach to this propagator at the eikonal approximation. In this case, the 
largest "+" component p + of the three-gluon vertex in the numerator shall cancel against 
a corresponding term in the denominator of the gluon propagator and the resulting Gf 
only contains its ?/ + -component. 

Similarly, we can draw the twist-6 amplitudes. There are many different TOPT di- 
agrams. Some typical diagrams of them are shown in Figs. 8 and 9 for the transverse 
and longitudinal coefficient functions, respectively. Figures 10 and 11 are corresponding 
geometric TOPT diagrams. It is really difficult to judge which diagrams give leading 
contributions to the nuclear size-dependent effect in Figs. 3, 4, 8 and 9. However, we 
shall show that the leading corrections in the corresponding geometric diagrams only 
come from a few of diagrams (Figs. 6a, 10a and 10b), which contain maximum number 
of colorless gluon pair, and such corrections can be simply resummed. 
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3 Normalization condition 

Any calculation of the DIS cross section needs a factorization scheme to separate the 
calculable partonic part from the nonperturbative matrix element. In such a scheme, 
the propagators with the spinor trace or the Lorentz indices linking two separated parts 
are decomposed. At the same time, an arbitrary constant may sneak into the undefined 
correlation function in this technique. Therefore, at first we determine the normalization 
constant in the definition of the correlation function. For this reason, we begin with the 
naive parton model, which defines the quark distribution functions f qi ( x B)- According to 
Fig. 12, the transverse structure function F T at the leading order takes the form, 

F T {x B ) = 2e^W^ 

= H e i / dx^Tr 

= ^2e 2 x B T 2q [x B ) = J2 e2 ifqA%B), (16) 

i i 

where we replace the propagator (k' 2 + ie)^ 1 with 5(k' 2 ) = x B 5(x — xb)/Q 2 due to the 
on-shell condition of the parton model. Thus, we have f qi (x) = 2xT 2q (x), where 

= J ^fe ixP+ y- < N\%(o)^My-)\N > . (17) 

The normalization constant in the gluon distribution function can be determined by 
using the contribution of gluons to F t (xb,Q 2 ) at the leading logarithmic approximation 
(Fig. 13). In the TOPT form, we have 

F^ LA (x B ,Q 2 ) = d^W^ 



xP +r y ■ nj^- 



Q 



2x B P+ 



7 • n7 



d^5(x - x B )^T 2q (x) 
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^ 2 f d 3 h 1 



A ^ p 
colour QTTQ-s 



E k s SE k s E P E ki \ E k s + E kl ~ E P 



x^7V[ 7 • A; l7 Q 7 • ks^dfTr^ ■ k 2l v l • M^T^s - — 



- E^/^|[(y) 2 + (l-(y) 2 )]/ 9 (y), 



where the momenta of partons in Fig. 13 are parameterized as 



y 'Q 2 



■2g 



(18) 



V 

k s 

h 



[yP + ,Q,yP" 



k 



{1 - x)P++ 2(l-x)P+ >- k± ' {1 - x)P " 



(19) 



The definition of T 2g (x) in Eq. (18) is 



T 2g (x) = J ^l^-e^v- < N\A +a (0)A^(y-)\N > d c A 

(20) 



or 



T 2g (x) = J ^f~e- p+ y- < N\F^(0)F^(y~)\N > (-^fdf, 

(21) 

where A +a — > F^n 13 / (xP + ) and F a/3 is the gluon field strength. The second definition 
has the gauge invariant form and presents singularity when x — ■> 0, which implies the 
exchange of soft gluon between the nonperturbative matrix and perturbative coefficient 
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function. However, all those soft gluons should be absorbed into the definition of parton 
distribution functions (PDFs) in the collinear factorization scheme. Therefore, the factor 
1/x in Eq. (21) should not give rise to any observable effects in the perturbation part of 
the amplitude. For avoiding misunderstanding, it is different from Ref. [6], we use the 
gluon distribution function 



P + dy 
2% 



< N\A +a {0)A +/3 {y~)\N> df . 



(22) 



12 



4 Geometric nuclear enhancement dynamics 

The nuclear target provides a special "laboratory" for the resummation of the higher 
twist effects. The contributions of the nuclear effects to the inclusive deep inelastic cross 
section are governed by the hadronic tensor. In the factorization scheme, the operator 
product in the hadronic tensor is expanded on the twist, where the higher twist ampli- 
tudes include the contributions of the multi-partonic configurations in the target. The 
contributions of the partons from different nucleons may give rise to a visible nuclear 
dependence in the structure functions. 

The nuclear effects in the structure functions originate from the correlations of partons 
among the bound nucleons. The parton correlation function is a generalization of the 
parton density beyond the leading twist. A general correlation function for two quarks 
and two gluons is defined as 



r P + dy P + dy 1 P + dy 2 ixiP + y - M x - Xl ) P + y - ^ x - X2)P + y - 
l2q29 ~ J 2tt 2tt 2tt 6 6 

X 2^I < Mn^)l + A +a ^)A + ^y^{y~)\A > df . (23) 

We always compare the structure function per bounded nucleon with that of a free 
nucleon to present the nuclear effects. Therefore, we do not expect a large nuclear depen- 
dence in the leading twist process, since 

T2 Q = J^^e- p+ y-^j<Amo) 1 ^(y-)\A> 

where \A > and \N > are the state vectors for nucleus and free nucleon, respectively. 
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In the general case, the exponential oscillations in the higher twist matrix element 
destroy any nuclear size enhancement that could come from the y~ integration [7,13]. 
The contributions of the integrations with respect to y{ and y^ to Eq. (23) are namely 
from the gluonic fields neighboring the targeted quarks. In this case, we assume that 



r P + dy P + dy 1 P + dy 2 iXlP+y ~ i(x _ Xl)P+y7 i(x _ X2)P+y - 1 
i 2i 2i 2t ' 2P+ 

xi < A|*(0-) 7 + A+%r)^feM^)l^ > df 

_ r P+dy- P+dyj P+dy? ^ P+y - ci(n)P+y7 i( x - aa )P+ y - 1 
J 2tt 2tt 2tt ' 2P+ 

x < NMO^+A+^y^A^iy^iy^lN > df . (25) 

Therefore, we do not expect a nuclear size effect in Eq. (25). However, some TOPT- 
diagrams for Ft(xb,Q 2 ) at twist-4 (Fig. 14) show that y7 = y^. We further consider 
the case where two gluonic operators consist a colorless cluster, which belongs to the 
same extra nucleon. In this case there is no strong correlation between the struck nucleon 
and the extra nucleon. The momentum exchange among the bounded nucleons can be 
neglected in comparison with a lager P + -scale. Thus, we assume that Xi ~ x 2 in Eq. (25). 
Now all nucleons listing on the integral path can contribute their gluons to the double 
scattering in Fig. 14 without distortion. Thus, we can write 



/ 



P+dy-P+dyi Jx , P+y -l 1 



2ir 2tt A 2P+ 

x < A\^(0-)^(y-)\A >< A\A +a (y^)A +(3 (y^)\A > df 
P+dy- P+d Vl JxiP+y - I ^ < ^(0-^+^-)]^ > 



2tt 2tt 2P+ 
x<A\A+«(y^)A+P(y^\A>d a / 
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AV3_i 




(26) 



R 2 



where 



P + dyj 
2ir 



<N\A +a (y^)A + ^y^\N>d' 



(27) 



The factors r/± and r)\\ are the transverse and longitudinal correlation functions, respec- 



tively. Q 2 rj± = Q 2 /R 2 can be schematically explained as the overlap probability of two 
partons, where 1/Q is the scale of a parton at momentum transfer Q 2 and R is the max- 
imum correlation length of two partons. Usually, R is regarded as an effective size of 
nucleon. The nuclear size effect occurs when the scattered quark exceeds the longitudinal 
size of a nucleon, and it reaches a maximum at xb = l/^itU^iv), where Ra is nuclear 
radius and the nucleon mass. 

Equation (26) shows a nuclear enhancement mechanism in the higher twist processes, 
i.e, the rescattering of the outgoing quark on the extra nucleons via the exchange of the 
contacting gluon-pair can give rise to a visible A 1//3 -dependence of the matrix element. 
We regard a pair of gluons, which are contacted in the y~-space as the effective particle 
(6 particle). 

From the geometric pictures in Figs. 6 and 7, a straightforward prediction is that the 
leading contributions to the nuclear size effect in the twist-4 structure functions are from 
Fig. 6a and its conjugate diagram. However, the oscillating factor can not be suppressed 
since Vi^y^-, although a quark pair in Fig. 6c can come from an extra nucleon and 
forms a colorless cluster. The contributions of Fig. 6b, 6d and 6e are also excluded since 
there is no contacting gluon pair, which can freely move along the longest path on the 
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y~-space. Thus, the leading contributions of the geometric effect at twist-4 are only from 
Figs. 6a and its conjugate diagram, with all the rest diagrams being sub- leading and 
negligible in the resummation. 

The leading diagrams in twist-6 amplitudes are Figs. 10a, 10b and their conjugates, 
which lead to the nuclear size effect ~ A 2 / 3 . On the other hand, the quark-gluon pairs 
in Figs. lOc-lOe are colored and they are strongly correlated with the other nucleons. 
Thus, we can not set X\ = x 2 in these diagrams to get the leading nuclear size effect. 
The integrated path of a contacting gluon-pair, which is connected with other nucleon 
via a gluon in Figs. lOf, is shorter than that in Fig. 10a. Therefore, we conclude that 
the contributions of Figs. lOc-lOf, lla-lld (and other diagrams, which have not dawn 
in Figs. 10 and 11) are < A 2 / 3 , and they are excluded in counting the leading effect. In 
consequence, we need resum only the leading contributions from such CLTOPT diagrams, 
where a maximum number of contracting gluon-pairs may freely move along in a larger 
nucleus. According to this rule, the contributions of the longitudinal structure functions 
are excluded in our resummation. Therefore, the geometry of the CLTOPT provides an 
ideal tool to choose the leading diagrams in the resummation of the geometric nuclear 
effect. 
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5 Geometric nuclear size effect in structure 

functions 

We calculate the contributions of higher twist amplitudes to the structure functions 
in this Section. We begin with the process shown in Fig. 14a. According to the results of 
the above section, this process provides a leading nuclear effect to the transverse structure 
function. Using the projection operator and TOPT rules, the contributions of Fig. 14a 
to the F T (x B , Q 2 ) at x\ — > x 2 are 



A a F T (x B ,Q 2 ) = 2e^W^ 

= ^47re-a: s J dx x dx 2 < - > colour d^ 



x-Tr 
4 



2(x 1 -x B )P+ Q 2 2x B P + 



xd°f5(x 2 - x B )^T 2q2 g 
= X>? / d*id*^ x ^*L-* B) M*uWA* - 1)9. (28) 

Similarly, we have the contributions of Fig. 14b to Ft at x 2 — > x\ 

A b F T (x B ,Q 2 ) = E^ 2 / d Xl dx 2 2 ™ s 2 Xb5{Xi ~ Xb) f qi ( Xl , Q*)(A^ _ 1)0 . (29) 
^ J Q l R l x 2 - x B 

Summing the contributions of Figs. 14a and 14b in the limit x\ ~ x 2 (see Eq. (A.l) in 
Appendix A) we get 



AF T (x B ,Q 2 ) =E^(^ 3 - l)Qx B ^-U(x B ,Q 2 ). (30) 

Here we point out that the results of hard part using the TOPT coincides with that using 
the covariant perturbation theory, which was given in [14]. 

17 



We emphasize that the contributions of the geometric nuclear effect to F^" lst ~ (xb, Q 2 ) 
are always negative at small xb due to the shape of the quark number density as a function 
xb- Equation (30) means that the double rescattering of an outgoing quark via an effective 
9 particle, which consists of a pair of contacting gluons, reduces the magnitudes of the 
transverse structure function. 

Now we consider the contributions of the higher twist amplitudes and sum over them. 
We calculate the processes in Fig. 15, where the cut line divides the initial gluons into two 
parts, each contains an even number of gluons. For simplicity, we indicate the particle 
as the dashed line. We have shown that these amplitudes provide maximum nuclear size 
effects in the structure functions at the twist-6 level. The contributions of Fig. 15a to 
Ft(xb, Q 2 ) at xi — > x 3 is 



A a F T<twist _ 6 (x B ,Q 2 ) = ^e 2 



dx\dx 2 dx 3 { 



2n:a s 2 x 2 B 5(x 3 - x B ) 



2x 3 T 2qig , (31) 



Q 2 R 2 (x 1 - x B ) (x 2 - x B ) 



where the matrix element is 




2x 3 T 2qig ( Xl ,Q 2 ) 
f P + dy~ P + dy 




(32) 



Therefore, we obtain 



A a F Tt twist-6(xB,Q ) 



Ee?/ d Xl dx 2 dx 3 (^) 2 (A^ - ifO 2 



x 2 B 5(x 3 - x B ) 



U(xi,Q 2 )- (33) 



(xi -x B )(x 2 -x B ) 
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Summing the contributions of three cut diagrams in Fig. 15 at the limit x\ — > xs, we get 
(see Eq. (A. 2) in Appendix A) 

AF T , twtst . 6 (x B ,Q 2 ) = E^(^) 2 (^ 1/3 " l^eV^/^.Q 2 ). (34) 

Similarly, we obtain the contributions from the twist- (2n + 2) amplitudes to F T (x B , Q 2 ), 
which contains the exchange of n-Q particles (see Eq. (A. 3) in Appendix A). The result 



is 

AF T , twlst ^ 2n+2) (x B ,Q 2 ) =1^(11)^/3 _ l) n ^ n x B £,f q AxB,Q 2 ). (35) 

Thus, we can easily write the sum of the contributions from all higher twist ampli- 
tudes, which give a leading nuclear size dependence of the transverse structure functions 
F T (x B ,Q 2 ), 

N 9-77-rv 1 r\ n 

AF T (x B ,Q 2 ) = T,e 2 U^nA^-lT^e n x B ^-U(x B ,Q 2 ) 



^Q 2 R 2 n\ "dx n B 

2na. , .Mr, _ x „ d 



cxp 



(A 1/s - 1)6x5— - 1 



f qi (x B ,Q 2 ). (36) 
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6 Geometric nuclear size effect in evolution 

equation 

In this section we discuss the contributions of the geometric higher twist effect to 
the DGLAP equation. We still work in the collinear infinite frame (see Eq. (2)). The 
coefficient functions with one quark loop corrections and at LLA(Q 2 ) can be modelled 
using either the standard scattering picture (Fig. 11) at p + ^> q~ or the dipole picture 
(Fig. 16) at q~ 3> p + (we call the latter frame as the dipole frame). There are different 
versions of the dipole model. A general dipole model for the DIS [15] is constructed in 
the /c^-factorization scheme, where the transverse momenta of the initial partons do not 
vanish and evolution dynamics is the BFKL equation [16]. In this work we shall show 
that the dipole picture can also be used in the collinear factorization scheme. 

As we have pointed out that the geometric nuclear effect is exposed in the collinear part 
of the outgoing quark interacting with the nuclear gluons, where the transverse momenta 
of the partons are negligible. Therefore, the geometric nuclear effect is hidden under the 
scattering picture Fig. 13, since the contributions of the transverse momentum flowing 
into the quark-gluon vertex in this picture is necessary in the derivation of the evolution 
(splitting) kernel. 

Different from the scattering picture, in the dipole picture (Fig. 16) the interacting 
part of quark with gluon (the dashed box in Fig. 16) can be separated from the splitting 
kernel using TOPT and then we may take the collinear approximation. At first step, we 
derive the DGLAP equation in the dipole model. The momenta of the partons in Fig. 16 
are parameterized as 
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Q2 _ 

q = q n - - — n, 
2q~ 



(37) 



i.e., 



E„ 



V2 



(q + + q ) ~ V2q = q' , 



(38) 



Thus, we have 



k 2 

k D ^ (z D q' + ± / _ ,k ± ,-z D q' ), 
izoq 

h ~ ((i - ^y- + fei -k ± , -(i - s y-)- 

A 1 — z D)q 



(39) 



On the other hand, 



p = (E p ,p ± ,p 3 ) = (y D P + ,0,y D P + ). 



(40) 



The relevant coefficient function which contains the splitting function at LLA(Q 2 ) in 
the dipole picture is 



/-rtwist—2 



d 3 h 



(27T)3 8 £r E k2 E q 



1 2 



X 



E k D - E k2 - E q\ 



k D 



Ee- fdk 2 ± dy D 1 



/ 



8vr 2 7 fc^ y D ^ D Q 2 



4 + (1 - z D ) 5 



twist— 2 



(41) 



where we introduce the probe function of the gluon 
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H 



twist— 2 




(42) 



and used the on-shell condition 



S((p + z D q) 2 ) 



1 



5(z D 



)■ 



(43) 



z D Q 2 



In the calculations of J2 \ M (7 — > gg)| 2 , we regard the virtual photon as an equivalent 
photon in the Weizsacker- Williams approximation [17] under condition q~ ^> p + . In this 
case, we consider the contributions of the transverse polarized photon and define zd- Note 
that the longitudinal coefficient function doesn't contribute to the DGLAP dynamics. 
The splitting function is defined as 



Using the factorized splitting function, we can construct a DGLAP equation in the dipole 
picture (see Appendix B and [18]) 




(44) 



i.e., 



P^ D ) = -[z 2 D + {l-z D )\ 



(45) 




(46) 



where x\ < yo < x B . 
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Comparing Eq. (46) with the DGLAP equation in the scattering picture, we have 
the following relations between the scattering and dipole pictures due to the different 
definitions about zd and zs~. 



P"(z s )^P q D Jz D ), 



(47) 



where zd = Ud/%b, %s — x b/us an d we use "D" and "S" to indicate the quantities in the 
dipole and scattering pictures, respectively. 

We emphasize that although both yo and ys are defined as the momentum fractions of 
the target carried by the gluon, they have different kinematical regions. It is interesting 
to note that a smaller value of xb corresponds to a smaller value of yo in the dipole 
model since xb = yo/zD', on the other hand, in the scattering picture the gluons can take 
a larger momentum fraction ys of proton-momentum even in the small Xb region since 
zs = xb/vs < 1- 

Using Eq. (47), one can rewrite Eq. (46) as a standard DGLAP equation 



where xb < ys < 1- 

In our dipole model the splitting kernel containing loop transverse momentum k± is 
separated from the probe function H twlst ~ 2 , in which we can take the collinear approxi- 
mation, i.e., neglecting the contributions of the loop transverse momentum. 

Now we discuss the corrections of the geometric nuclear effect to the DGLAP equation 
using the dipole picture. The nuclear gluons can collinearly attach to two quark legs in 




(48) 
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any ways. However, the leading contributions to the geometric nuclear effect are from 
following probe functions: 



H t W ist-A = H t^st-2 27Tas{A l/3 _ l)_^|_e ^Lf g (y D ,Q*), 



y D Q 2 R 2 dy D 



(49) 



at twist-4 (see Fig. 17), and 



H twist~6 = H twist - 2 (27ca s ) 2 (A l/3 - 1 



2 X B 



1 d 



®^^f 9 (yD,Q 2 ). 



y 2 D Q*R* 2! dyl 



(50) 



at twist-6 (see Fig. 18). 

We sum the leading contributions of the nuclear size-effect at all order of twist to the 
DGLAP equation, the result is 

df qi (x B ,Q 2 ) 
dlnQ 2 

_ a s f 1 dy s s f ^u(„ n2] 



+- 



« s f XB dy D D y D 



L 



2tt J x% y D mK x B 



)\x B 



y D =VD 



(51) 



We don't need to consider the nuclear size-effect on the gluon legs since these diagrams 
are excluded as the sub-leading contributions in our resummation. 
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7 Discussions and summary 

One of our results, Eq. (36), is consistent with the Qiu-Vitev resummation for 
AF T (x B ,Q 2 ) in Ref. [6], only differing from an unknown constant. The reason is that 
these two approaches are based on different dynamics. The derivation of the Qiu-Vitev 
formula assumes that the leading contributions to the nuclear size effect in the structure 
functions are arisen from the singularities of the gluon correlation functions. This argu- 
ment does not indicate a lot of neglecting graphs, which give the sub-leading contributions. 
In this respect, the geometry of the CLTOPT improves the Qiu-Vitev derivation. 

It is interesting to compare Eq. (51) with the Mueller-Glauber model [2]. These 
two models are written in the dipole picture but using different factorization schemes. 
Besides, the Glauber approach of the multiple scattering theory is based on the eikonal 
approximation, which assumes that the scattering on different nucleons are incoherent. 
Thus, the phase shifts from each scattering can sum up to an exponential form. In the 
TOPT-language, the eikonal approximation implies that all propagators connecting with 
nuclear gluons take their forward component. On the other hand, Eq. (51) provides a 
different picture, where the backward propagator alternates with the forward one when 
scattered quark passes through the nucleus. 

It is well known that there is an obvious nuclear screening effect in the structure 
function F 2 A at small xb [19]. As we have pointed out that the nuclear size effect in the 
longitudinal stricture function at leading order vanishes. The nuclear screening effect at 
small xb region can be expressed in the ratios of the structure functions as 

2 , F Ai2 (x B ,Q 2 ) 



R(x B ,Q 2 ) = 



F N , 2 (x B ,Q 2 ) 
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x B F A)L (x B , Q 2 ) + x B F AtT (x B , Q 2 
F N , 2 (x B ,Q 2 ) 

~ + F N:2 (x B , Q 2 ) ■ (52) 
From these results, one can fix the parameter G. 

The nuclear PDFs are one of fundamental knowledge for understanding the experi- 
mental data at the Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider 
(LHC). The PDFs are calculable if the nonperturbative input distributions and their QCD 
evolution dynamics are known. Therefore, the corrections of nuclear circumstances to the 
evolution equations and input distributions become an active topic. The nuclear PDFs 
can be modelled using the DGLAP equation from the relating DIS data [20]. The results 
show that the input PDFs in the bound nucleon differ by the shadowing factors from 
that in the free nucleon. However, this result does not means that the corrections from 
the higher twist effects in such researches are unimportant since the shadowing factors 
are phenomenological. In fact, the nonlinear modification of the gluon recombination to 
the DGLAP equation as a higher twist effect can't be neglected even in DIS on the free 
nucleon [21]. As mentioned already in works [6], the higher twist resummation may be 
of importance in the extraction of the nuclear input PDFs. Our results remind us of 
that, except the contributions of the gluon recombination effect to the DGLAP equation, 
the geometric higher twist effect in the nuclear input PDFs and QCD evolution equation 
should be considered in extracting the nuclear parton distributions from the measured 
DIS cross sections. 

As we have mentioned that the higher twist effects and their summation may provide a 
valuable information for the multi-parton correlations at long distance in the target. The 
geometric nuclear effect presents an interesting freedom in high energy nuclear physics— 
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the contacting gluon-pair, which is soft at the hadronic level since it carries almost zero- 
momentum (xi — x 2 — 0), however, the effect is really calculable in the perturbative 
QCD. 

In summary, the higher twist effect in deep inelastic scattering on the nuclear target 
are studied using time ordered perturbation theory. We have shown that the collinear 
rescattering of the outgoing quark on the extra nucleon via the contacting gluon-pair 
dominates a nuclear size dependent effect. The leading contributions of nuclear-enhanced 
effect to the structure function and DGLAP evolution equation are resummed by using 
the geometric properties of the higher twist amplitudes. 

Acknowledgments: This work was supported by National Natural Science Foundations 
of China 10075020, 50193013, and 10475028. 
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Appendix A: 

We list some algebraic identities in the derivations in Section 5. 
In F T) twist-A 



when A — > 0. 

In FTjwist-6 



dxidx 2 
— / dxi 



5(x 2 - x B ) 5(x 1 - x B ) 



Xi - x B x 2 — x b 
5(xi — x b + A) 5(x 1 -x b ) 



S(x 2 - x 1 - A) 



-J dx 1 [S'(x 1 -x B )\ 



dxt 



(Al) 



dx±dx 2 dx 3 



5(x 3 - x B ) 



+ 



5(x 2 - x B ) 



+ 



8(xi - x B ) 



_(xi - x B )(x 2 - x B ) (x± - x b )(x 3 - x b ) (x 2 - x B )(x 3 - x B ) 
x5(x 2 — xi — A)5(x 3 — x 2 — A) 



dx\dx 2 dx 3 



5(x 3 -x B ) 8(x 2 -x B ) , 5(x 1 -x B ) 



2A 2 



2A 2 



-/ 

= - j dx 1 dx 2 — [S'(x 1 - x B ) - S'(x 2 - x B )]5(x 2 - x x - A) 
= - J dx 1 — [8'(x 1 - x B ) - 5\xi -x B + A)] 
= J dx-^8"(xi - x B ) 



5(x 2 — xi — A)5(x 3 - x 2 



1 d 2 
2\dx 2 B 



when A — > 0. 

In Fxjwists 



(A2) 
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/ 



dx\dx 2 dx 3 dx 4 



5(x A - x B ) 



+ 



8(x 3 - x B ) 



(xi - X B ) (X2 - x B ) (^ 3 - x B ) (xi - x B )(x 2 - x B )(x 4 - X B ) 
5(x 2 -x B ) 5(x 1 -x b ) 



(xi - x B ){x 3 - x B )(x 4 - x B ) (x 2 - x B ){x 3 - x B )(x A - X B ) 
xS(x 2 — x\ — A)b[x 3 — x 2 — A)5(x 4 — x 3 — A) 



-/ 



dx\dx 2 dx 3 dx 4 



S(x 4 - x B ) S(x 3 - xb) _ S(x 2 - xb) 8(xi - x B ) 



2A 3 



6A 3 



6A 3 2A 3 
xS(x 2 — x\ — A)5(x 3 — x 2 — A)5(x 4 — x 3 — A) 

= - J dx 1 dx 2 dx 3 -^^[5' \x 3 - x B ) - 5'(x 2 - x B ) - S'(x 2 - x B ) + 5'(xi - x B )} 

x5{x 2 - Xi- A)5{x 3 - x 2 - A) 

= - / dx 1 — [5"(x 1 - x B + A) - S"( Xl - x B )} 

= - J dxi^5"'(xi - x B ) 



1 d 3 
3!cte|' 



(A3) 



when A — > 0. 
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Appendix B: 

In this Appendix we derive a set of complete DGLAP evolution equations in the dipole 
frame. At first step, we give the definitions of the parton distributions in this frame. The 
QCD evolution equations of the PDFs are irrelevant to the probe, therefor, we may use a 
virtual gluon or a virtual quark, which connects with an extra color source "S" , to replace 
the virtual photon in DIS. Thus, the first order of the DIS amplitudes in the dipole frame, 
can be expressed as the four graphs shown in Fig. 19. 

We take Fig. 19a as an example. The cross section of a virtual gluon scattering off a 
proton in the dipole frame and in the collinear factorization scheme is 



For convenience, we neglect the symbol "D" for the dipole picture. According to the par- 
ton model, 



where C a is a function describing the probe vertex. In the above equation and below, we 
will use the momenta of the partons to represent the particles entering cross sections and 
matrices. The gluon distribution in the proton at the first order QCD can now be defined 

as 





da a (qp -> fci) = C a 5(y - x B ), 



(B.2) 



fgi.X B ) 



1 

a 



a 



da a (qP -> hX). 



(B.3) 



Similarly, we have three more definitions corresponding to Fig. 19b-19d. 
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fgM = ^rda b (qP - hX), (BA) 



b 



f q (x B ) = ^rda c (qP -> fciX), (5.5) 



and 



/,(zb) = ^rd<j d (qP -> fcxX). (5.6) 
w 

On the other hand, we can write the formulae for the four cross sections of Fig. 19 in 
the TOPT form as 



5 p 



|M,(P->pX:)| : 



Ep — Er, — E 



v - 



2E n 



-r-r d 3 k X 

y (2tt) 3 25 x 



1 H 3 ki 
^ M a , b (p - gfc! - p)(27r) 4 5 4 (p - g - h)- " 



SEpE q 



(27T)32E fcl ' 



(5.7) 



and 



da c , d (qP -> fciX) 



|i|M 9 (P^X)| 2 



i2 r 



5p — E p — Ex 



2E„ 



n 



X- 



1 



8E p E q 



M c , d (p - gfcx - p)(27r) 4 <5 4 (p - g - fcx) 



(2vr) 3 25 x 



(5.8) 



(27r) 3 25 fcl - 

Comparing Eqs. (B.7) and (B.8) with (B.l), we get the definitions of the quark and gluon 
distributions 
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„, VJ „.. M -^\M g (P^pX)\ 2 

Hip 



f g {y)dy 



i 2 r 



Ep — E p — Ex 



2E„ 



-pr d 3 k x 
1 x 1 (2tt) 3 2£; x ' 



(5.9) 



/ g (^ = fW(P-pX)| 2 



i?P — E„ — Ex 



2E n 



yr d 3 k X 

^^ir^Ex' 



(5.10) 



and the probe vertices 



d<jj(qp — > fci) 



— MAp 

Oj OHipHiq 



qki — > p)(27r) 5 (p — q — ki 



5(y-x B ) 



d 3 ki 



(2nf2E kl 



(5-11) 



in the TOPT form, where j — a,b,c and d. 

When the scale Q 2 is increased by a small amount due to the radiation of a parton, 
i.e., Q 2 — > Q 2 + AQ 2 , the corresponding change of the PDFs, for example in Fig. 20a' is 



dfg(x B ) 



1 



C, 



Hip 



pX)f 



Ep — E p — Ex 



2E„ 



^jS^E- x Ha{qP 



qp) 



= J dyf q (y)H a (qp -> qp), 



(S.12) 



where 



5 a (<?P -> gp) 
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a 



■da a (p -> z D qh -> p)- E r\M(q -> fcfc 2 



5,: — ^fa ~~ 



25r. 



d 3 A; 2 



= -<5(^ - -2-) 



y . dfc^ a. 



P 9 9 ( Z d), 



(2n) 3 2E k2 

(5.13) 



y x B fci 27r 

where P gg (zr)) is the splitting function for g — > gg. The momenta of the partons in Fig. 
20 are parameterized as Eqs. (39) and (40). From Eqs. (B.12) and (B.13), we get a 
DGLAP evolution equation 



and 



df g (x B ,Q 2 
d\nQ 2 



a s f XB dy 2 y 



P gg {z D ) = 2C 2 (G) 



1 - z D z D 

+ V z D (l - z D ) 



(5.14) 



(5.15) 



1 - Z£> 

Along similar lines, we derived the evolution equations corresponding to Figs. 20b' 



20d' 



df Qi {xB,Q 2 ) _ <**_ f XB dy f ( 2x p , 
dlnQ 2 "2tt4 y J9W^ qi gK XB h 



(5.16) 



and 



dlnQ 2 



df g {x B ,Q 2 ) _ a^X [ x *dy , 2 ,y_, 
- o / 2 /<z;UW ) r a<n\ )■> 

Z7T ■ 1 Jxf, y X_b 



(5.17) 



df qi (x B ,Q 2 ) 
dlnQ 2 

The splitting functions are 



= 2^4 7 4(2/ ' g)P ^y 



(5.18) 



= g [4 + (1 - zdY 



(5.19) 
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1 + (1 - z D ) 
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P m (zD)=C 2 (R) 



(5.20) 



and 



P qiqj (z D ) = C 2 (R) 



1 - Z D 



(5.21) 



Note that two partons indicated k in Fig. 20 are on mass-shell in the TOPT framework, 
therefore, one can extract the probe function }j twlst ~ 2 from the amplitude (Fig. 20b') and 
generalize it to the higher twist processes as in Section 6. 

The probabilistic form of the splitting functions in Eqs. (B.15) and (B.21) have 1/(1 — 
z D ) singularities, which are caused by the emission of a low-momentum gluon. Physically, 
the infrared singularities can be cancelled when the virtual diagrams, for example Fig. 21 
are added to the real contributions. As we have pointed out in Ref. [22], an advantage of 
TOPT is that it exposes the relation between the real and virtual diagrams without an 
actual calculation of the virtual diagrams. In fact, from the TOPT-decomposition of the 
diagrams, we can easily get 



only half of the probe-vertex connects with the partonic matrix in the virtual diagrams. 




That is 



34 



The contributions from others virtual diagrams can be derived by similar method. In 
consequence, we get the evolution equations 



dU(*B,Q 2 ) 
dlnQ 2 



2n Jxl 



x b dy 



f g (y,Q 2 )P qtS 



x B 



(5.24) 



and 



df g (x B ,Q 2 ) 
dlnQ 2 
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Figure Captions 
Fig. 1 Elemental vertices of QCD at the collinear limit. 
Fig. 2 Elemental vertices of QED at the collinear limit. 

Fig. 3 The twist-4 TOPT diagrams at lowest order for F t (xb, Q 2 ), the conjugant diagrams 
have not drawn. 

Fig. 4 As similar to Fig. 3, but for F L (x B ,Q 2 )- 

Fig. 5 An example of time-space picture corresponding to Fig. 3a. 

Fig. 6 Space-time pictures of the amplitudes corresponding to Fig. 3. 

Fig. 7 Space-time pictures of the amplitudes corresponding to Fig. 4. 

Fig. 8 A part of twist-6 TOPT diagrams at lowest order for Ft(xb,Q 2 ), the conjugant 

diagrams have not drawn. 

Fig. 9 A part of twist-6 TOPT diagrams at lowest order for Fl(xb,Q 2 ), the conjugant 
diagrams have not drawn. 

Fig. 10 Space-time pictures of the amplitudes corresponding to Fig. 8. 

Fig. 11 Space-time pictures of the amplitudes corresponding to Fig. 9. 

Fig. 12 The naive parton model, which defines the quark distributions. 

Fig. 13 The leading order contribution of gluon to F t (xb,Q 2 ) in the standard parton 

(scattering) model. 

Fig. 14 The leading perturbative contribution to the nuclear size-dependent of F t (xb, Q 2 ) 
at twist-4, where an outgoing quark rescattering on the extra nucleon via the contacting 
colorless gluon-pair exchange. 

Fig. 15 The leading diagrams at the twist-6 level, where using the dashed line to denote 
a contacting colorless gluon pair. 
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Fig. 16 Schematic diagrams for the decompositions of a DGLAP equation at the dipole 
frame in TOPT. 

Fig. 17 Space-time picture for the probing function at twist-4, which gives the leading 
contribution to the nuclear size-dependent of the DGLAP equation. 
Fig. 18 Similar to Fig. 17, but at twist-6. 

Fig. 19 The definitions of the parton distributions in the dipole frame. Where S is the 
color source. 

Fig. 20 The one-loop corrections to the parton distributions in the dipole picture. 
Fig. 21 The virtual diagrams corresponding to Fig. 20a. 
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Fig.9 
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Fig.13 
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